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ABSTRACT 


This  report  describes  the  continued  evaluation  of  the  19- 
element  Alaskan  Long  Period  Array  (ALP A),  which  was  conducted  by  Texas 
Instruments  Incorporated,  at  the  Seismic  Data  Analysis  Center  over  the  per¬ 
iod  1  April  1972  to  31  March  1973. 

The  major  areas  of  study  in  the  evaluation  were: 

•  Noise  analysis 

•  Regionalization  of  seismic  events 

•  Matched  filter  performance 

•  Analysis  of  S-wave  processing  for  the  Kurile  Islands  - 

Kamchatka  region 

•  Seismic  event  detection  thresholds 

•  Behavior  of  seismic  discriminants 

A  total  of  379  events  were  processed  and  analyzed  in  the  course 
of  this  evaluation.  Where  applicable,  earlier  ALPA  data  and  results  are  dis¬ 
cussed  in  conjunction  with  the  present  results. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  necessarily  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical 
Applications  Center,  or  the  US  Government. 
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SECTION  I 

|  INTRODUCTION 

This  report  presents  the  results  of  a  continuation  of  the 
evaluation  of  the  19-element  Alaskan  Long  Period  Array  (ALPA).  It  extends 
the  analysis  reported  in  Special  Report  No.  4,  Extended  Array  Evaluation 
Program  (Heiting,  et  al,  1972).  Emphasis  was  placed  on  building  a  large  data 
base  so  that  detection  methods  and  discrimination  parameters  of  seismic 
events  might  be  studied  on  a  regional  basis.  Specific  areas  of  investigation 
include: 

•  Noise  analysis 

•  Regionalization  of  seismic  events 

•  Matched  filter  performance 

•  Analysis  of  S-wave  processing  for  the  Kurile  Islands  - 
Kamchatka  region 

•  Seismic  event  detection  threshholds 

•  Behavior  of  seismic  discriminants 

Data  from  previous  years  has  been  included  in  the  data  base  and,  where 
applicable,  is  used  in  the  evaluation. 

The  data  base  and  the  data  processing  methods  are  described 
in  Section  II.  Details  of  the  analysis  of  events  from  specific  areas  of  interest 
are  discussed  in  Sections  III  through  VIII.  Section  IX  summarizes  results, 
presents  conclusions  and  suggests  possible  areas  of  further  analysis  utilizing 
the  ALPA  array. 

Mr 
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SECTION  II 


DATA  BASE  AND  ANALYSIS  METHODS 

The  results  presented  in  this  report  are  based  on  a  compilation 
of  seismic  events  and  presumed  explosions  that  were  recorded  in  1970,  1971, 
and  1972.  The  event  parameters  are  listed  in  Appendix  A.  Each  event  is 
named  by  a  three  part  designator  consisting  of  a  three  letter  abbreviation 
for  the  region,  the  Julian  date,  and  the  hour  (GMT)  of  occurrence.  These 
three  parts  are  separated  by  symbols  indicating  the  year  of  occurrence: 
slashes  denote  1970  events,  asterisks  denote  1971  events,  and  dashes  denote 
197  2  events. 

A.  EVENT  SELECTION 

The  method  of  selecting  events  to  be  processed  from  those 
recorded  in  1972  was  to  compile  a  list  of  all  events  having  epicenters  in  the 
general  area  of  interest  from  the  available  event  lists.  These  were:  the 
Preliminary  Determination  of  Epicenters  Monthly  Summary  (PDE) ,  the 
SDAC/LASA  Weekly  Summary  (LASA),  the  NORSAR  Seismic  Event  Summary 
(NORSAR),  and,  for  the  period  20  February  to  19  March  1972,  the  Internation¬ 
al  Seismological  Month  list  (ISM)  provided  by  Massachusetts  Institute  of  Tech- 
nology  (Lincoln  Laboratories).  When  more  than  one  of  these  lists  reported  a 
given  event,  preference  of  choice  of  epicentral  data  was  in  the  order  ISM, 

PDE,  LASA,  NORSAR.  Appendix  B  breaks  down  the  disposition  of  the  events 
proposed  for  processing. 

Two  major  event  suites  were  processed  for  1972.  The  first,  a 
winter  suite,  was  composed  of  events  recorded  during  the  period  1  January  1972 


through  20  March  1972.  The  second,  a  summer  suite,  was  composed  of 
events  recorded  during  the  period  1  June  1972  through  31  August  1972.  This 
separation  in  time  of  the  two  suites  was  chosen  to  permit  an  investigation  of 
possible  differences  in  detection  capability  of  ALPA  in  summer  and  winter. 

Other  events  outside  this  time  period  were  processed,  either  to  use  as 
reference  waveform  matched  filters  or  to  build  up  the  presumed  explosion 
data  base. 

The  data  base  for  this  report,  including  1970  and  1971  events, 
totaled  524  seismic  events  and  32  presumed  explosions.  It  was  felt  that  the 
1970  partial-array  events  could  be  included  in  this  data  base,  since  full-array 
and  partial-array  surface-wave  beamforming  gains  have  an  average  difference  of 
only  one  dB  (Heiting,  et  al,  1972).  Therefore,  there  should  be  no  appreciable 
difference  in  detection  or  measured  surface-wave  magnitudes  between  full- 
array  and  partial-array  events. 

The  data  base  is,  therefore,  composed  of  the  following  time 

periodsr 

Number  of  Number  of 

Seismic  Events  Presumed  Explosions 


1970 

61 

7 

1971 

96 

13 

Winter,  1972 

165 

3 

Summer,  1972 

20  2 

9 

The  histogram  (Figure  II- 1)  shows  the  number  of  events  from 

each  information  source  and  the  total  number  processed  as  a  function  of  m  . 

b 

Note  that  the  PDE  events  predominate  at  higher  values  of  m  ,  and  LASA  and 

b 

NORSAR  events  at  lower  values  of  m^.  The  ISM  events  are  fairly  evenly  dis¬ 
tributed  throughout  the  range  of  values  of  m  . 
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DISTRIBUTION  BY  DATA  SOURCE  OF  EVENTS  PROCESSED 


B.  DATA  PROCESSING  METHODS 
1.  Signal  Processing 

The  array  configuration  of  ALPA  is  shown  in  Figure  II-2.  in 
ali  processing  of  seismic  signals  discussed  in  this  report,  the  raw  data 
recorded  by  the  triaxial  seismometer  at  each  array  were  rotated  by  means  of 
a  transformation  of  coordinates  to  form  three  orthogonal  components  of  ground 
motion,  one  vertical  and  two  horizontal.  The  two  horizontal  components  are 
oriented  such  that  the  radial  component  always  lies  in  the  direction  of  the 
great  circle  path  of  the  event  in  question  and  the  transverse  component  always 
lies  perpendicular.  Therefore,  the  Love  wave  energy  will  always  occur  on  the 
transverse  (T)  component  and  the  Rayleigh  wave  energy  on  the  vertical  (V)  and 
radial  (R)  components. 

A  beamsteered  trace  was  then  formed  for  each  component  of 
motion  (T,  V,  and  R)  using  all  good  sites.  The  number  of  good  sites  for  each 
event  is  given  in  Appendix  A.  the  velocities  used  in  beamforming  were  4.  0 
km/sec  for  the  Love  wave,  3.  5  km/sec  for  the  Rayleigh  wave,  and  distance- 
dependent  velocities  for  the  shear  waves. 

The  beamsteered  traces  were  filtered  using  the  standard 
0.025-0.055  Hz  bandpass  filter,  an  appropriate  reference  waveform  matched 
filter,  and  five  chirp  filters.  A  second  bandpass  filter  of  0.  033  to  0.083  Hz 
was  applied  to  the  July  and  August  events  to  investigate  the  frequency  depend¬ 
ence  of  the  surface-wave  magnitude  measurements. 

During  the  period  6  June  1972  through  28  August  1972,  new  fil¬ 
ter  amplifiers  were  installed  at  all  19  sites  of  ALPA,  changing  the  quantiza¬ 
tion  level  from  0.  565  millimicrons  per  computer  count  (m/i/cc)  to  0.  28 
m^/cc.  Table  II- 1  lists  the  dates  when  the  work  was  begun  and  completed  at 


TABLE  II- 1 

AMPLIFIER  CHANGE-OVER  TIMES  FOR  ALPA 


SITE 

DATE  BEGUN 

1 

29  June 

2 

20  July 

3 

06  June 

4 

22  August 

5 

24  June 

6 

1 3  July 

7 

19  August 

8 

19  August 

9 

19  August 

10 

28  August 

11 

23  August 

12 

23  August 

13 

23  August 

14 

17  August 

15 

17  August 

16 

17  August 

17 

2  0  July 

18 

20  July 

19 

12  July 

DATE  COMPLETED 

06  July 

20  July 
06  July 

22  August 
06  July 

18  July 

19  August 
19  August 
19  August 
28  August 

23  August 
23  August 
2  3  August 
18  August 

17  August 

18  August 

21  July 
28  August 
13  July 


II- 6 


each  site.  For  any  given  clay  during  this  period,  some  sites  were  operating 
at  the  old  (0.  565  m/i/cc)  quantizing  level  and  some  at  the  new  (0.  28  mfx/cc). 
This  was  corrected  for  in  the  edit  process,  where  those  sites  operating  at 
0.  28  m /x/cc  were  scaled  to  the  0.  565  m/Lt/cc  quantization  level. 

2.  Noise  Data  Base 

The  raw  data  forming  a  noise  sample,  as  recorded  by  the 
triaxial  seismometers,  were  rotated  to  form  three  components  of  ground 
motion  vertical,  north-south,  and  east-west.  The  three  components  of  the 
single-site  noise  data  were  plotted  for  three  or  four  sites  of  the  array  to  deter¬ 
mine  whether  there  were  any  seismic  signals  present  which  were  not  reported 
by  any  of  the  information  sources.  If  a  signal  was  found,  a  new  time  period 
was  sought  for  the  noise  sample.  These  plots  also  allowed  a  check  for  spikes 
in  the  noise  sample.  Sites  which  had  anomalously  high  or  low  power  levels 
(as  computed  by  the  edit  routine)  were  dropped  from  the  analysis.  The 
remaining  good  sites  were  then  used  to  compute  cross -power  spectral  ma¬ 
trices  from  which  the  average  RMS  noise  levels  could  be  measured.  The 
resulting  noise  data  base  is  listed  in  Appendix  C.  Finally,  frequency- 
wavenumber  spectra  were  computed  to  investigate  the  source  azimuths  of 
the  peak  micro  seismic  noise. 
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SECTION  III 
NOISE  ANALYSIS 


A.  INTRODUCTION 

This  analysis  of  the  ambient  noise  field  at  ALPA  is  an  exten¬ 
sion  of  the  analysis  performed  for  the  preceeding  year  (Heiting,  et  al,  1972). 
The  objectives  of  this  analysis  are  to  characterize  the  noise  field  in  terms 
of  the  RMS  noise  level  and  the  directionality  of  the  noise  and  to  determine 
whether  the  results  are  consistent  with  those  obtained  during  the  previous 
year. 

One-hour  noise  samples  were  taken  at  approximately  ten-day 
intervals  throughout  the  year,  as  listed  in  Appendix  C.  All  data  were  re¬ 
sampled  to  a  two-second  sampling  interval  and  divided  into  256-second  (128- 
point)  segments.  The  data  were  examined  for  sites  and  segments  which  were 
dead  or  contained  spikes  or  glitches.  These  bad  sites  and  segments  were 
dropped  from  further  analysis.  Next,  a  crosspower  matrix  was  generated 
for  each  noise  sample  at  52  frequencies  from  0.  0  to  0.  2  Hz  (Af  =  .  00391  Hz) 
by: 

•  Removing  the  mean  from  the  data  of  each  component  at  each 
site 

•  Discrete  Fourier  transforming  the  individual  data  segments 

•  Hanning  the  transforms 

•  Cross -multiplying  to  obtain  the  crosspower  terms 

•  Stacking  over  all  segments 


In  order  to  make  comparisons  with  the  results  of  the  previous 
year,  the  spectra  were  not  corrected  for  instrument  response.  The  nominal 
value  of  0.  565  millimicrons  per  computer  count  (m/i/cc)  at  25  seconds  (0.  04 
Hz)  was  used  to  normalize  the  power  density  spectra.  Samples  taken  after 
24  June  1972  were  scaled  appropriately  to  keep  the  normalization  value  at 
0.  565  m/z/cc,  as  described  in  Section  II. 

B.  RMS  NOISE  LEVELS 

Figures  III  - 1 ,  III-.',  and  III- 3  illustrate  the  behavior  of  the 
RMS  noise  level  for  the  vertical,  east-west,  and  north-south  components 
respectively,  as  observed  at  ALPA  during  the  period  May  1971  through  Dec¬ 
ember  1972.  (The  1971  data  was  taken  from  Special  Report  No.  4,  Extended 
Array  Evaluation  Program.  )  These  RMS  noise  levels  were  computed  from 
the  average  across  the  sites  over  the  one-hour  time  gate.  The  bandwidth  18 
to  40  seconds  (0.  025  to  0.  055  Hz)  was  used  in  calculating  the  RMS  noise 
levels. 

Inspection  of  the  data  shows  the  RMS  noise  levels  of  the  three 
components  follow  the  same  general  pattern.  There  are,  however,  a  few  ex¬ 
ceptions.  For  example,  day  190  of  1971  shows  a  high  RMS  noise  value  on  the 
north- south  component  with  much  lower  values  on  the  vertical  and  east-west 
components.  The  power  spectra  for  this  noise  sample  indicate  that  this  is 
probably  due  to  non-propagating  long  period  noise,  which  is  higher  on  the 
north- south  component  than  on  the  other  two  components. 

The  high  RMS  noise  values  have  been  previously  attributed  to 
coherent  noise  with  a  spectral  peak  at  about  .  055  Hz  or  long -period  non¬ 
propagating  noise  with  a  fairly  flat  spectrum  in  the  range  0.  0  to  0.  04  Hz.  An 
examination  of  the  1971  and  1972  noise  data  indicate  that  the  high  RMS  noise 
levels  from  day  240  to  day  320  of  1971  and  at  days  60,  220,  and  350  of  1972 
are  due  to  long-period  non-propagating  noise.  Furthermore,  the  high  noise 
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levels  at  days  330  and  361  of  1972  are  at  least  partly  due  to  long-period  non¬ 
propagating  noise.  However,  the  high  noise  levels  at  day  322  of  1971  and 
days  1  and  320  of  1972,  as  evidenced  by  the  power  spectra,  do  not  contain 
long -period  non-propagating  noise  and  are  due  to  coherent  microseismic 
noise.  This  coherent  noise  may  be  storm-generated. 

An  interesting  feature  of  the  noise  data  is  found  in  the  RMS 
noise  level  data  for  the  period  day  241  of  1972  through  day  361  of  1972.  These 
data  indicate  an  upward  trend.  This  trend  is  in  contrast  to  the  corresponding 
period  in  1971,  where  the  data  imply  a  fairly  constant  background  RMS  noise 
level  punctuated  by  short-duration  bursts  of  higher -energy  noise.  To  deter¬ 
mine  whether  this  is  a  real  difference  in  the  behavior  of  the  noise  field  would 
require  more  noise  samples  from  the  day  241  to  day  361  period  and  an  ex¬ 
tension  of  this  period  to  approximately  day  100  of  1973. 

The  following  conclusions  may  be  drawn  from  the  RMS  noise 

level  data: 

•  With  the  exception  of  the  unexplained  upward  trend  in  the  data 
at  the  end  of  1972,  the  RMS  noise  values  appear  to  have  a  fair¬ 
ly  constant  background  level  throughout  the  year,  ranging  be¬ 
tween  7  and  10  millimicrons. 

•  The  higher  RMS  noise  levels  for  the  most  part  appear  to  be 
due  to  bursts  of  long-period  non-propagating  noise  superim¬ 
posed  on  the  background  noise  level.  A  few  of  the  high  RMS 
noise  levels  may  be  due  to  storm  generated  noise. 

C.  DIRECTIONALITY  OF  THE  NOISE 

The  source  azimuths  of  microseismic  noise  were  measured 
on  high-resolution  frequency-wavenumber  spectra.  Each  of  the  frequency- 
wavenumber  spectra  was  computed  at  the  frequency  corresponding  to  the 
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maximum  value  of  the  power  density  spectrum  for  that  sample.  This  peak 
frequency,  which  ranged  from  0.  047  to  0.  063  Hz;,  was  used  because  it  is 
generally  highly  coherent.  Azimuths  were  computed  only  from  vertical- 
component  data,  since  the  horizontal  components  contain  both  Love  and  Ray¬ 
leigh  energy.  The  peak  frequencies,  azimuths,  and  velocities  are  listed  in 
Appendix  C. 

The  source  azimuths  of  microseismic  noise  as  recorded  at 
ALPA  are  shown  in  Figure  III-4.  In  this  figure,  the  source  azimuth  of  the 
peak  microseismic  noise  is  indicated  by  an  arrowhead.  The  source  azimuths 
of  the  clearly  discernable  secondary  peaks  (2  dB  or  less  below  the  primary 
peak)  are  indicated  by  circles.  The  range  of  source  azimuths  for  the  con¬ 
tinuum  of  energy  6  dB  or  less  below  the  primary  peak  is  indicated  by  a  line. 

The  results  of  this  analysis  agree  with  the  results  of  the  ana¬ 
lysis  of  the  previous  year  (Heiting,  et  al,  1972).  The  predominant  micro¬ 
seismic  noise  source  azimuth  during  1972  lies  in  the  range  of  azimuths  125° 
to  150°,  which  coincides  with  the  western  Canada  and  United  States  coastlines. 
A  secondary  range  of  source  azimuths  was  found  between  180°  and  230°,  which 
coincides  with  the  Cook  Inlet  region  of  Southern  Alaska.  A  possible  explana¬ 
tion  for  these  apparent  microseismic  noise  sources  is  that  storm-generated 
waves  are  channeled  into  the  Alexander  Archipelago  along  the  western  coast 
and  the  Cook  Inlet,  resulting  in  the  release  of  relatively  large  amounts  of 
wave  energy  along  restricted  coastal  areas. 

In  summary,  the  results  of  the  analysis  of  the  ambient  noise 
field  at  ALPA  agree  with  those  of  the  preceeding  year.  High  levels  of  long- 
period  noise  occur  only  in  the  winter  months.  The  source  azimuths  of  micro¬ 
seismic  noise  as  recorded  at  ALPA  rarely  coincide  with  azimuths  to  the  area 
of  interest. 
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SECTION  IV 

REGIONALIZATION  OF  SEISMIC  EVENTS 


The  data  base  was  regionalized  by  mapping  the  epicenters  of 
all  available  events.  Each  set  of  closely  groupe  i  events  was  then  considered 
to  form  the  population  of  a  seismic  region.  This  decomposition  of  the  data 
base  into  separate  regional  populations  made  possible  an  investigation  of 
matched  filter  responses,  detection  thresholds,  and  standard  discriminants 
in  more  detail  than  was  previously  possible. 

The  regions  so  chosen  are: 

•  Southern  Kurile  Islands  (KUR)  This  region  contains  events 

with  epicenters  in  or  near  the  Kurile  trench  having  epi- 
central  distances  (A)  to  ALPA  of  approximately  35°-40°. 

•  Kamchatka  (KAM)  This  region  contains  events  with  epicenters 

along  the  eastern  coast  of  Kamchatka,  near  the  Komandor- 
sky  Islands,  and  in  the  northern  end  of  the  Kurile  trench. 
These  events  have  epicentral  distances  of  about  25°-30°. 

•  Central  Asia  (CA)  This  region  contains  events  with  epicenters 

in  Sinkiang,  Tibet,  and  the  Hindu  Kush.  These  events 
have  epicentral  distances  of  about  60°-80°. 

•  Caspian  Sea  (CSP)  This  region  contains  events  with  epicenters 

near  the  Caspian  Sea.  It  includes  events  from  the  Cauca¬ 
sus  Mountains  and  northwestern  Iran.  These  events  have 
epicentral  distances  of  about  70°-80°. 
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Southern  Iran  (SIR)  This  region  contains  events  with  epicenters 
near  the  coast  of  southeastern  Iran.  These  events  have 
epicentral  distances  of  about  85°. 

Greece- Western  Turkey  (GT)  This  region  contains  events  with 
epicenters  in  Greece,  Italy,  the  Adriatic  Sea,  the  Aegean 
Sea,  and  Western  Turkey.  These  events  have  epicentral 
distances  of  about  70°-80°. 

Eastern  Kazakh  Test  Area  (EKZ)  This  region  contains  events 
which  have  epicenters  in  the  Eastern  Kazakh  test  area  and 
which  are  all  presumed  explosions.  These  events  have 
epicentral  distances  of  approximately  60°. 

Taiwan  (TWN)  This  region  contains  events  which  have  epi¬ 
centers  in  or  near  Taiwan  and  the  southern  Ryukyu  Islands. 
These  events  have  epicentral  distances  of  approximately 
70°. 

The  locations  of  these  regions  are  shown  in  Figure  IV- 1.  The 
luded  in  each  of  these  regions  are  so  noted  in  Appendix  A. 


MILLER  MODIFIED  MERCATOR  PROJECTION 
MAP  SCALE*  0.750  IN./  15  DEG.  LONGITUDE 

FIGURE  IV -1 


SECTION  V 

MATCHED  FILTER  PERFORMANCE 


A.  INTRODUCTION 

Matched  filters  were  applied  to  long-period  signals  as  recorded 
at  ALPA  to  evaluate  their  effectiveness  in  increasing  the  signai-to -noise  ratio 
of  dispersed  seismic  signals.  Both  the  reference  waveform  matched  filters 
and  chirp  filters  were  applied  to  the  transverse  Love  wave  and  the  vertical  and 
radial  Rayleigh  waves  of  a  test  event.  The  goals  of  this  analysis  were: 

•  To  determine  potential  signal-to-noise  ratio  gains  of  reference 
waveform  and  chirp  filters 

•  To  compare  the  relative  effectiveness  of  reference  waveform 
and  chirp  filters 

•  To  evaluate  the  effectiveness  of  matched  filters  in  increasing 
the  surface-wave  detection  capability  of  ALPA. 

Matched  filter  performance  was  analyzed  in  terms  of  signal-to- 
noise  ratio  improvement  over  the  equivalent  bandpass  signal-to-noise  ratio. 
Each  signal-to-noise  ratio  was  calculated  as  the  ratio  of  the  peak  value  of  the 
signal  waveform  to  the  RMS  value  of  the  noise  measured  in  a  gate  ahead  of  the 
signal.  The  signal-to-noise  ratio  improvement  of  a  matched  filtered  beam  over 
the  corresponding  bandpass  beam,  expressed  in  decibels,  is: 

Improvement  (dB)  =  20  log  (S^matchecj/SNRBp) 

Since  the  signals  are  not  noise -free,  the  signal  amplitudes  are  actually  signal 
plus  noise  amplitudes.  For  this  reason,  we  will  refer  to  the  signal  plus  noise- 
to-noise  ratio  (SNNR)  from  this  point  on. 
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B„  REFERENCE  WAVEFORM  MATCHED  FILTER  RESULTS 

A  suite  of  27  reference  waveform  events  was  selected  for  this 
evaluation  of  reference  waveform  matched  filtering.  The  approximate  loca¬ 
tions  of  these  events  are  shown  in  Figure  V-l.  The  name  and  associated 
parameters  for  the  event  corresponding  to  each  numbered  location  in  the  fig¬ 
ure  are  given  in  Table  V-l.  The  selection  criteria  were:  good  SNNR,  shallow 
focus  (less  than  60  km),  and  location  in  an  area  where  no  reference  waveform 
event  has  been  previously  selected.  The  length  of  the  reference  waveforms 
was  chosen  in  the  following  manner:  for  events  at  large  epicentral  distances, 
the  length  was  selected  to  include  possible  multipath  energy,  since  small 
changes  in  event  epicenter  location  would  not  be  expected  to  significantly 
change  the  multipath  structure.  This  situation  is  reversed  for  events  at 
small  epicentral  distances;  for  such  events,  small  changes  in  event  epicenter 
location  could  significantly  change  the  multipath  structure.  Therefore,  the 
lengths  of  reference  waveforms  having  small  epicentral  distances  were  cho¬ 
sen  so  as  to  exclude  any  possible  multipath  energy. 

The  SNNR  improvements  obtained  from  reference  waveform 
matched  filtering  of  77  events  from  1972  are  given  in  Table  V-2,  along  with 
the  reference  waveform-test  event  separation.  (There  is  no  duplication  of 
test  events  for  different  reference  waveform  filters.  )  The  letters  in  paran- 

thesis  to  the  left  of  the  name  of  each  reference  waveform  refer  to  the  region 
(Section  IV)  to  which  it  belongs. 

Considering  those  reference  waveforms  for  which  there  are 
enough  test  events  to  make  discussion  meaningful,  the  following  statements 
can  be  made  about  the  behavior  of  individual  reference  waveform  filters. 

•  Reference  waveform  KAM/242/00  yielded  good  Love  wave  im¬ 

provement  and  poor  Rayleigh  wave  improvement.  The  mean 
Rayleigh  wave  improvements  are  not  as  good  as  those  obtain¬ 
ed  last  year  using  a  different  set  of  test  events. 
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TABLE  V-l 

REFERENCE  WAVEFORM  EVENT  DATA 


Event 

Name 

KOM-171  -01 
KAM*2  06*03 
KAM/  242  /  00 
KUR-063-23 
KUR/ 219/01 
SAK*251*16 
K  YU/  2  06/22 
SIB/156/10 
ERS-  165-10 
ERS/ 241/  14 
PIP-039-03 
CHI/212/13 
BUR/210/10 
USM-057  -2 3 
KAZ/249/04 
SIN-002 -10 
SIN-047-23 
SIN-084-08 
SIN*  2 1 9* 1 5 
TAD-077 -09 
HIN-  176-15 
IRA/242/ 16 
IRA-101-02 
PAK-028  - 1 0 
GRC/184/00 
ADB/231/02 
CRS/287/05 


Date  Time  Lat. 

M/ D/ Y  (hr-min-sec)  (°N) 

06/19/72  01-43-48  54.4 

07/31/71  03-45-05  52.6 

08/30/70  00-38-40  52.1 

02/03/72  23-10-40  50.2 

08/07/70  01  -43-19  43.8 

09/06/71  16-59-53  48.0 

07/25/70  22-41-11  32.2 

06/05/70  10-31-54  63.4 

06/13/72  10-45-05  54.9 

08/29/70  14-59-23  51.1 

02/08/72  03-37-52  19.3 

07/31/70  13-10-47  28.6 

07/29/70  10-16-19  26.0 

02/26/72  23-31  -09  50.6 

09/06/70  04-02-57  49.8 

01/02/72  10-27-35  42.8 

02/16/72  23-19-20  41.7 

03/24/72  08-1  1  -53  42.9 

08/05/71  15-21  -53  36.  1 

03  17/72  09-17-10  40. 1 

06/24/72  15-29-22  36.2 

08/30/70  16-17-31  37.4 

04/10/72  02-07-28  33.2 

01/28/72  10-26-54  26.6 

07/03/70  00-41-01  38.7 

08/19/70  02-01-53  41.4 

10/14/70  05-59-57  73.3 


Lon. 

(°E) 

Depth 

(km) 

mb 

Location 

number 
Fig.  V-l 

168.  6 

33 

5.  0 

1 

160.  7 

33 

4.  5 

2 

i  j 6 

33 

5.  2 

3 

155. 7 

33 

4.  5 

4 

148.  3 

33 

5.  0 

5 

143.  0 

16 

5.  9 

6 

131.  7 

34 

6.  1 

7 

146.  2 

33 

5.  5 

8 

126.  4 

33 

4.  8 

9 

135.  3 

33 

5.  4 

10 

122.  0 

50 

5.  7 

11 

103.  6 

25 

5.  5 

12 

95.  4 

59 

6.  5 

13 

97.  3 

33 

5.  3 

14 

78.  1 
87.  3 
80.  7 
87.  4 
77.  7 
69.  7 
69.  7 
56.  0 
56.  6 
66.  3 
20.  4 
19.  8 
55.  1 
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REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENTS 

(PAGE  1  OF  5) 


1 

dB  SNNR  Improvement  Over 

Equivalent  Bandpass  Filter 

Event 

RMF/  Test 

1  (.  025  -  .  055 

Hz) 

Designation 

Separation  (km) 

T 

V 

R 

(KAM)  KUR-063-23 

(RMF) 

KUR-2  35  -  14  AL 

55 

9.  1 

7.  7 

4.  7 

KUR-194-00QC 

89 

5.  0 

4.  2 

4.  3 

KUR-2  18-22AL 

132 

3.  4 

4.  5 

1.  6 

KUR-2  09-00AL 

221 

3.  2 

10.  3 

8.  3 

KUR  -154-01  AL 

286 

1.  9 

4.  0 

4.  3 

KUR-2  1  6-02QC 

426 

-0.  5 

2.  3 

2.4 

Mean  Improvement  for  KUR-063-23 

3.  7 

5.  5 

4.  3 

Standard  Deviation 

3.  2 

2.  9 

2.  3 

(KAM)  KOM-  171-01 

-4-  (RMF) 

KOM-  183-02AL 

175 

0.  3 

3.  0 

3.  1 

KOM-229  - 1 OQC 

210 

2.  1 

1.  5 

4.  5 

KOM- 1 8  0  -  04  AL 

289 

2.  8 

3.  6 

3.  6 

KOM  -  153-21TD 

303 

5.  9 

3.  6 

1.  5 

KAM  - 1 86  -  1  3  AL 

366 

3.  8 

1.  3 

3.  3 

KAM-22  9-2  1  AL 

397 

3.  6 

1.  4 

2.  9 

Mean  Improvement  for  KOM-171-01 

3.  1 

2.  4 

3.  1 

Standard  Deviation 

1.  9 

1.  1 

1.  0 

(KAM)  KAM*2 06*03 

4- (RMF) 

KAM  - 1 7  3  -00QC 

49 

-1.  3 

-5.  3 

-1.7 

KAM  - 1 80  - 14  AL 

49 

4.  0 

0.  6 

1.  0 

KAM  - 1  92  -  12  AL 

130 

2.  6 

1.  9 

1.  0 

KAM- 17  3-1  OQC 

157 

1.7 

-3.  7 

-0.  7 

KAM  - 17  7  -  1 7  AL 

162 

-1.  0 

1.  6 

0.  5 

K  AM-17  9-06  AL 

192 

3.  1 

-2.  6 

-2.  1 

KAM-193-08QC 

306 

-1.  3 

-2.  7 

-2.  8 

KAM  - 1  68  -09AL 

378 

1.  8 

3.7 

-1.0 

KAM-1  58-1 0AL 

418 

-1.4 

1.  9 

-1.8 

KAM  - 1  57  -04  AL 

429 

1.  8 

1.  0 

2.  6 

Mean  Improvement  for  KAM*206*03 

1.  0 

-0.  6 

-0.  5 

Standard  Deviation 

2.  0 

2.  7 

1.  6 

V-5 
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REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENTS 
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Event 

Designation 

RMFy  Test 
Separation  (km) 

dB  SNNR  Improvement  Over 
Equivalent  Bandpass  Filter 
(.  025  -  .  055  Hz ) 

T 

V 

R 

(KAM)  KAM/  242/  00 

>(RMF) 

■■ 

KAM-23  1  -  1 9  AL 

102 

5.  5 

0.  2 

KAM-156-07  AL 

147 

4.  0 

-1.7 

KAM-233-08QC 

164 

4.  2 

-3.  5 

KAM- 199- 08  A  L 

322 

6.  0 

0.  6 

w 

Mean  Improvement  for  KAM/242/00 

4.  9 

m 

0.  3 

Standard  Deviation 

1.  0 

2.  1 

(SIR-CSP)  IRA- 101  -( 

32  4-(RMF) 

IRA- 187  - 1  6  AL 

497 

-0.  7 

6.  8 

IRA-216-22  AL 

557 

-0.  2 

-2.  3 

MM 

IRA- 155  -08  AL 

579 

-0.  2 

0.  8 

-0.  1 

IRA-196-13QC 

648 

5.  6 

6.  2 

7.  8 

IRA-184-12  AL 

648 

2.  0 

5.  5 

3.  7 

IRA-182-17QC 

667 

0.  6 

1.  3 

IIQ-175-08AL 

969 

-1.8 

1.  6 

IIQ-166-04  AL 

977 

4.  0 

-4.  3 

Mean  Improvement  for  IRA- 101 -02 

1.  2 

1.4 

1.  9 

Standard  Deviation 

2.4 

4.  3 

3.  9 

(TWN)  PIP-039-03 

(RMF) 

TWN-212-16AL 

223 

2.  5 

7.  2 

7.  0 

T  WN-  178-08QC 

267 

-0.  2 

10.  1 

9.  8 

TWN-160-09AL 

274 

-1.0 

10.  9 

13.  6 

TWN- 1 98-13  AL 

494 

-0.  6 

6.  6 

4.  1 

TWN-182-18TD 

563 

0.  6 

3.  4 

4.  7 

RYU-155-02  AL 

591 

0.  0 

3.  8 

3.  3 

RYU-J  97-02  AL 

631 

0.  5 

8.  8 

10.  5 

Mean  Improvement  for  PIP-039-03 

5.  8 

7.  6 

Standard  Deviation 

1.  1 

3.  2 

3.  8 

V-6 


TABLE  V -2 

REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENTS 

(PAGE  3  OF  5) 


Event 

Designation 

RMF/ Test 
Separation  (km) 

(EKZ)  KAZ/249/04 

-<-(RMF) 

EKZ-345-04AL 

0 

EKZ-159-01QD 

7 

EKZ  -307-01  AL 

78 

Mean  Improvement  for  KAZ/249/04 
Standard  Deviation 

PAK-028-1  0 

-4-  (RMF) 

PAK- 162-1 1  AL 

179 

PAK-179-06QC 

522 

PAK- 179-1  0AL 

522 

IRA-22 1  -1  9  AL 

550 

Mean  Improvement  for  PAK-028-10 
Standard  Deviation 

CHI/212/13 

-4-  (RMF) 

TSI-  156-23  AL 

797 

TIB-198-02  AL 

854 

TSI-243  - 1 5  AL 

1117 

■  ■  -  t 

- — - — - i 

dB  SNNR  Improvement  Over 
Equivalent  Bandpass  Filter 
(.  025  -  .  055  Hz) 


T  V  R 


(CA)  BUR/210/10 

BUR-179-09AL 
TSI-  180-03  AL 


■  (RMF) 


(CA)  SIN-047 -2 3 

SIN-154-06AL 
SIN-187-01  AL 


(RMF) 
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TABLE  V-2 


REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENTS 

(PAGE  4  OF  5) 


Event 

Designation 

RMF/  Test 
Separation  (km) 

dB  SNNR  Improvement  Over 
Equivalent  Bandpass  Filter 
(.  025  -  .  055  Hz) 

T 

V 

R 

(CA)  SIN-084-08 

SIN- 187  -04  AL 

SIN- 192-1  9AL 

4-  (RMF) 

88 

112 

1.4 
-0.  7 

5.  4 

1.6 

1.9 

1.  1 

(CA)  HIN-176-15 

HIN-178-20AL 
AFG-2  11-17  AL 

4-  (RMF) 

67 

492 

3.  6 

1.8 

0.  7 

4.  5 

-1.2 

3.  3 

(CA)  SIN-002-10 

TSI-154-16AL 

4-  (RMF) 

856 

-1.6 

-3.  1 

0.4 

(CA)  TAD-077-09 

AUB-181-03QC 

-4- (RMF) 

198 

2.  0 

-2.4 

1.  5 

(CSP)  IRA/242/ 16 

IRA-202  -  1  3  AL 

IIQ- 164-1  3  AL 
IIQ-165-00AL 

-4-  (RMF) 

179 

1000 

1000 

5.  1 
-2.4 

0.  6 

6.  7 

0.  5 
-1.  1 

6.4 

4.  3 

3.4 

(GT)  GRC/ 184/00 

MED-  187-18AL 
TUR-2 1 6-21  AL 

-4-  (RMF) 

222 

1079 

0.  8 

2.  7 

0.  7 

5.  6 

1.  5 

5.4 

(GT)  ALB/231/02 

YUG-243-00AL 

-4- (RMF) 

412 

2.4 

3.  0 

0.  8 

(KUR)  KUR/  219/01 

KUR-232-00AL 

4-  (RMF) 

34 

2.  0 

4.  8 

3.4 

CRS/287/06 

NVZ -241-05  AL 

4- (RMF) 

0 

2.  5 

10.  6 

6.  5 

- -  - - _ 
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TABLE  V-2 

REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENTS 

(PAGE  5  OF  5) 


Event 

_ 

Designation 

RMF/  Test 
Separation  (km) 

dB  SNNR  Improvement  Over 
Equivalent  Bandpass  Filter 
(.  025  -  .  055  Hz) 

T 

V 

R 

SAK*251*16 

KUR-2 1 1  -2  1 AL 

(RMF) 

978 

1.6 

0.  5 

0.  6 

ERS/ 241/14 

ERS-165  -1 OQD 

'♦■(RMF) 

729 

1.9 

4.  2 

0.  9 

SIB/156/10 

SIB-238-04  AL 

-4- (RMF) 

912 

0.4 

3.  2 

3.  9 

ERS-  165-10 

ERS-222-2  0AL 

-4- (RMF) 

217 

2.  5 

1.  o 

1.  o 

USM-057  -23 

CRS-244  -14  AL 

-4- (RMF) 

230 

5.  6 

2.  9 

3.  0 

KYU/  2  06/22 

RYU-209-  16AL 

-4-  (RMF) 

764 

1.  1 

6.  0 

5.  8 

I 

I 

I 

1 

1 
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Reference  waveform  KUR-063-23  yielded  good  Rayleigh  wave 
improvements  for  all  test  events  and  good  Love  wave  improve¬ 
ments  on  all  but  the  last  two  test  events.  This  decrease  in 
Love  wave  improvement  appears  to  be  due  to  increases  in  ref¬ 
erence  waveform-test  event  separation. 

Reference  waveform  KOM-171-Ol  yielded  fairly  good  Rayleigh 
and  Love  wave  improvements.  The  average  improvements  for 
the  three  components  are  essentially  the  same. 

Reference  waveform  KAM*206*03  yielded  poor  Rayleigh  and 
Love  wave  improvements.  This  indicates  that  this  reference 
waveform  is  not  representative  of  the  events  occurring  near 
its  epicentr al  location.  Therefore,  this  reference  waveform 
should  be  replaced. 

Reference  waveform  IRA-101-02  yielded  poor  Rayleigh  and  Love 
wave  improvements  on  five  test  events  and  good  Rayleigh  wave 
improvements  on  three  test  events  (located  near  each  other). 
Therefore,  there  should  be  another  reference  waveform  avail¬ 
able  in  this  area  to  apply  to  the  group  of  events  which  show 
poor  improvements  when  IRA- 101 -02  is  applied. 

Reference  waveform  KAZ/249/04  yielded  good  Rayleigh  wave 
improvements  and  poor  Love  wave  improvements.  The  Ray¬ 
leigh  wave  improvements  are  better  than  those  of  the  previous 
year. 

Reference  waveform  PAK-028-10  yielded  good  Rayleigh  wave 
improvements.  Love  wave  improvements  were  poor  except 
for  the  first  test  event,  which  is  much  closer  to  the  reference 
waveform  than  the  other  three  test  events. 


•  Reference  waveform  PIP-039-03  yielded  poor  Love  wave  but 

excellent  Rayleigh  wave  improvements  with  the  exception  of 
the  last  test  event. 

Thus,  with  the  exception  of  reference  waveform  KAM*Z06*03, 
reference  waveform  matched  filtering  yields  good  SNNR  improvements  (>  3 
dB)  on  at  least  one  surface  wave  propagation  mode. 

The  effect  of  reference  waveform-test  event  separation  upon 
SNNR  improvement  is  illustrated  by  Figure  V-2.  A  straight  line  least-mean- 
square-error  fit  to  the  data  points  and  the  95  per  cent  confidence  limits  for 
this  fit  are  shown.  The  plot  indicates  that  the  SNNR  improvement  obtained 
by  the  reference  waveform  matched  filters  decreases  gradually  with  increas¬ 
ing  reference  waveform-test  event  separation.  The  fitted  line  has  a  slope  of 
-3  dB  per  1000  km  of  separation. 

A  measure  of  the  variation  in  SNNR  improvement  among  test 
events  for  a  given  reference  waveform  matched  filter  is  the  standard  devia¬ 
tion  of  the  improvement,  which  is  also  listed  in  Table  V-2  under  the  mean 
SNNR  improvement  values  for  each  reference  waveform.  The  standard  de¬ 
viations  of  the  event  groups  are  generally  rather  large,  having  values  com¬ 
parable  to,  or  in  some  cases  greater  than,  the  mean  values.  This  indicates 
that,  for  most  of  the  event  groups,  the  idea  of  an  "average"  improvement  may 
be  meaningless.  Since  the  primary  reason  for  computing  an  average  im¬ 
provement  is  to  obtain  a  correction  factor  for  the  apparent  surface  wave  mag¬ 
nitude  of  an  event  detected  by  a  matched  filter,  an  event  group  with  a  large 
standard  deviation  in  improvement  may  produce  an  erroneous  magnitude  es¬ 
timate.  For  example,  the  two-standard-deviation  uncertainty  in  Love  wave 
magnitude  for  Kamchatka  events  matched  filtered  by  KUR-063-23  is  +_  0.  3 
magnitude  units. 

The  effectiveness  of  reference  waveform  matched  filters  was 
estimated  (using  1972  LR-V  data)  by  computing  the  percentage  of  SNNR 
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improvements  less  than  0  dB,  between  0  and  2  dB,  etc.  These  percentages 
are:  20%  of  the  improvements  were  less  than  0  dB,  30%  were  between  0  and 
2  dB,  20%  were  between  2  and  4  dB,  15%  were  between  4  and  6  dB,  9%  were 
between  6  and  8  dB,  and  6%  were  greater  than  8  dB.  Therefore,  even  though 
the  standard  deviations  are  so  large  as  to  render  "average"  SNNR  improve¬ 
ment  values  almost  meaningless,  it  is  possible  to  say  that  the  probability 
that  the  output  trace  of  a  reference  waveform  matched  filter  will  show  a  pos¬ 
itive  SNNR  improvement  is  0.  80  and  the  probability  that  the  improvement  will 
be  between  2  and  8  dB  is  approximately  0.  45.  Thus,  for  approximately  half 
of  the  events  to  which  reference  waveform  matched  filters  are  applied,  the 
SNNR  improvement  will  be  large  enough  to  aid  detection  studies. 


C.  CHIRP  FILTER  RESULTS 

Linear  chirp  matched  filters  were  applied  to  the  beamsteered 
signal  outputs  of  76  of  the  events  which  were  used  in  the  reference  waveform 
matched  filter  evaluation.  These  chirp  filters  were  specified  and  applied  in 
the  frequency  domain,  using  a  chirp  bandpass  of  0.  025  to  0.  055  Hz.  After 
application  of  these  chirp  filters,  the  data  were  inverse-transformed  to  ob¬ 
tain  time  domain  chirp  filter  outputs. 

The  chirp  filter  response  function  is: 


G{K)  =  e 


i2  7r(C/  N)  (K  -  K  )‘ 


if  K  <  K<  K 
Li  H 


where: 


G(  -K)  =  G(K)* 


if  0<K<  or  Kh<  K<N/2 


K  is  the  discrete  Fourier  transform  frequency  index 

KL  and  KH  are  the  lowest  and  highest  frequencies  in  the  passband 
Kq  is  the  frequency  index  at  which  zero  phase  shift  occurs 
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N  is  the  number  of  transform  points 

C  is  a  parameter  which  controls  the  length  of  the  correspond¬ 
ing  time  domain  waveform. 

This  yields  a  dispersive  time  domain  waveform  with  a  linear 
group  delay  and  essentially  flat  amplitude  at  all  periods  in  the  band  corres¬ 
ponding  to  Kl<K<Kh  (Harley,  1971). 

Five  chirps  were  applied  to  each  test  event.  Their  lengths  were 
centered  about  the  assumed  optimum  length  and  differed  by  increments  of  ±  50 
seconds.  The  assumed  optimum  length  was,  in  each  case,  picked  from  the 
plots  of  chirp  length  vs.  distance  as  presented  in  the  report  of  the  preceding 
year  (Heiting,  et  al,  1972).  The  SNNR  improvement  for  the  test  event  was 
then  measured  from  the  best,  in  terms  of  amplitude  and  shape,  of  the  five 
chirp  responses. 

Table  V-3  presents  the  chirp  filter  SNNR  improvement  in  dB 
over  the  equivalent  bandpass  filter  of  these  76  events.  The  mean  SNNR  im¬ 
provements  and  corresponding  standard  deviations  are  also  listed  for  each 
region. 

Considering  those  regions  which  contained  four  or  more  test 
events  (three  in  the  case  of  EKZ ),  the  following  comments  on  chirp  matched 
filter  results  can  be  made. 

•  KAM  -  SNNR  improvements  were  fair  on  all  three  components. 

The  Rayleigh  wave  improvements  are  slightly  better  than  the 

Love  wave  improvements. 

•  CA  and  CSP  -  SNNR  improvements  were  poor  on  Love  waves 

and  fair  on  Rayleigh  waves. 

•  SIR  -  SNNR  improvements  were  fair  on  Love  waves  and  good 

on  Rayleigh  waves. 
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TABLE  V- 3 


CHIRP  FILTER  IMPROVEMENTS 
(PAGE  1  OF  3) 


Event 

Designation 

Region 

dB  SNNR  lmproveme 
Equivalent  Bandpass 
(.  025  -  .  055 

;nt  Over 

Filter 

Hz) 

T 

V 

R 

KOM-  153-21TD 

RAM 

-0.  9 

-0.  3 

4.  1 

KUR-1  54-01  AL 

RAM 

0.  7 

0.  8 

1.  3 

KAM-156-07AL 

RAM 

-0.  1 

2.  3 

-0.  3 

RAM- 157-  04  AL 

RAM 

2.  1 

3.  9 

2.4 

RAM- 158-  10AL 

RAM 

0.  0 

7.  8 

4.  7 

RAM-  165  -  04  AL 

RAM 

0.  2 

1.  0 

-0.4 

RAM- 168- 09  AL 

RAM 

3.  6 

3.  7 

2.  3 

RAM  -  173-OOQC 

RAM 

0.  6 

-0.  1 

0.  5 

RAM- 173  -  1  OQC 

RAM 

2.  2 

1.  8 

2.  6 

RAM- 1 77-17  AL 

RAM 

2.  6 

2.  5 

0.  3 

RAM- 179 -06  AL 

RAM 

0.  9 

3.  5 

2.  3 

ROM  - 1 8  0  -  04  AL 

RAM 

2.  0 

5.  5 

6.  5 

RAM-  180-  14  AL 

RAM 

5.2 

0.  0 

2.  1 

ROM-  183-02  AL 

RAM 

1.2 

4.  0 

6.4 

RAM- 186-  1  3AL 

RAM 

4.  6 

-0.  2 

2.  7 

RAM-192-12AL 

RAM 

-0.  2 

-1.5 

0.  0 

RAM- 193- 08  QC 

RAM 

2.  9 

1.  9 

1.  9 

RUR-194-00QC 

RAM 

4.4 

6.4 

5.  2 

RAM  -  1 99  -08  AL 

RAM 

1.4 

1.  8 

1.6 

RUR-209-00AL 

RAM 

0.  6 

3.  0 

3.4 

RUR-2 1 6-02QC 

RAM 

1.  1 

3.  7 

1.8 

RUR-218-22AL 

RAM 

4.  0 

6.  0 

1.7 

ROM- 22  9-  1  OQC 

RAM 

3.  0 

3.4 

3.  7 

RAM-229-21  AL 

RAM 

1.  7 

0.  7 

2.  1 

RAM-231  - 1 9  AL 

RAM 

4.  7 

2.  5 

3.  2 

RAM-23  3- 08  AL 

RAM 

2.  7 

0.  3 

0.  1 

RUR-235  -  14  AL 

RAM 

5.  7 

4.7 

4.  8 

Mean  Improvements  for  Ramchatka 

2.  1 

2.  6 

2.  5 

Standard  Deviation 

1.  8 

2.  2 

1.  9 

RUR-232  -  2  3  AL 

RUR 

1.  1 

5.  5 

3.4 

SIN-154-06AL 

CA 

2.  3 

4.  2 

5.  0 

TSI-154-16AL 

CA 

-0.  3 

0.4 

2.  0 

TIB  -  170-  04  AL 

CA 

-1.  3 

-0.  7 

1.  0 
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TABLE  V- 3 

CHIRP  FILTER  IMPROVEMENTS 
(PAGE  2  OF  3) 


Event 

Designation 

Region 

dB  SNNR  Improvement  Over 
Equivalent  Bandpass  Filter 
(.  025  -  .  055  Hz) 

T 

V 

R 

HIN-  178-20AL 

CA 

2.  3 

4.  2 

1.6 

TSI-180-03AL 

CA 

2.4 

1.4 

3.  1 

AUB-  181  -03QC 

CA 

2.  5 

-2.  4 

-1.  9 

SIN-187-01AL 

CA 

-1.1 

4.  2 

4.  1 

SIN-187-04  AL 

CA 

2.  8 

0.  2 

-0.  1 

SIN-192-19AL 

CA 

-1.9 

2.  8 

1.  5 

KTB-206- 14  AL 

CA 

3.  8 

0.  6 

-0.  8 

AFG  -211-17  AL 

CA 

3.  8 

2.  1 

4.  0 

Mean  Improvements  for  Central  Asia 

1.4 

1.  5 

1.8 

Standard  Deviation 

2.  1 

2.  2 

2.  2 

IIQ-164-13AL 

CSP 

-1.  3 

1.6 

1.  5 

IIQ-165-00AL 

CSP 

-1.8 

1.  1 

1.  3 

IIQ-166-04  AL 

CSP 

0.  2 

-0.  1 

0.  1 

IIQ-17  5-08  AL 

CSP 

1.9 

1.2 

2.  3 

IRA-202  -  1  3  AL 

CSP 

0.  0 

4.4 

4.  2 

Mean  Improvements  for  Caspian  Sea 

-0.  2 

mm 

wgm 

Standard  Deviation 

1.4 

■SHI 

EH 

IRA- 155  -08  AL 

SIR 

1.4 

-0.  5 

0.  1 

IRA-182-17QC 

SIR 

2.  8 

4.  7 

4.7 

IRA- 184- 12  AL 

SIR 

0.  0 

4.  2 

3.  3 

IRA- 187-16  AL 

SIR 

-1.  1 

4.  9 

7.  3 

IRA-196-13QC 

SIR 

3.  0 

3.  6 

5.  0 

IRA-216-22  AL 

SIR 

4.  5 

5.  7 

4.  5 

Mean  Improvements  for  Southern  Iran 

1.8 

3.  8 

4.  2 

Standard  Deviation 

2.  1 

2.2 

2.4 

MED-  187-18AL 

GT 

2.  2 

2.  0 

2.  5 

TUR-2  16-21  AL 

GT 

0.  5 

2.  6 

2.  3 

YUG-243-00AL 

GT 

0.  4 

2.  4 

0.  2 

RYU- 1 55  -02  AL 

jjjggpf  amg 

2.4 

m 

1.  3 

TWN-160-09AL 

4.  5 

3.  5 

TWN-  178-08QC 

HHBExSHU 

4.  2 

6.  1 

6.  2 

TABLE  V-3 

CHIRP  FILTER  IMPROVEMENTS 
(PAGE  3  OF  3) 


dB  SNNR  Improvement  Over 

Equivalent  Bandpass  Filter 

Event 

_ 

025  -  .  055 

Hz) 

Designation 

Region 

T 

V 

R 

TWN-  182-18TD 

TWN 

2.  1 

6.  7 

5.  6 

RYU- 1 97-02  AL 

TWN 

2.  9 

4.  3 

5.  3 

TWN-198-.13AL 

TWN 

3.  4 

3.  8 

2.  5 

T  WN-2  12  - 1 6  AL 

TWN 

1.  0 

4.  8 

5.  5 

Mean  Improvement  for  Taiwan  Region 

2.  9 

4.  4 

4.  3 

Standard  Deviatior 

1.  2 

1.  7 

1.8 

EKZ  - 1 59-01QD 

Eastern  Kazakh* 

4.  6 

4.  1 

6.  1 

EKZ  -307-01  AL 

Eastern  Kazakh’1' 

0.  6 

5.  4 

5.  5 

EKZ  -  345  -04  AL 

Eastern  Kazakh* 

2.  5 

6.  8 

6.  5 

Mean  Improvement  for  Eastern  Kazkah 

2.  6 

5.  4 

6.  0 

Standard  Deviation 

2.  0 

1.4 

0.  5 

NVZ  -241  -05  AL 

Novaya  Zemlya’ 

4.  5 

8.  5 

7.  1 

TSI-156-23AL 

Tsinghai  Prov. 

2.  0 

0.  0 

0.  8 

ERS-165-10QD 

Eastern  Russia 

0.  6 

5.  3 

0.  9 

TIB  -  1  98  -02  AL 

Tibet 

-2.  0 

-0.  3 

- - 

-2.  1 

RYU-209-16AL 

Ryukyu  Islands 

3.  5 

7.  3 

6.  9 

KUR-2  11-21  AL 

Kurile  Islands 

2.  9 

2.  0 

1.8 

IRA-22 1  - 19 AL 

Eastern  Iran 

-3.  1 

3.  9 

3.  7 

ERS-222  -  20  AL 

Eastern  Russia 

2.  0 

4.  3 

6.  0 

SIB-238-04AL 

Siberia 

1.  9 

4.  5 

3.  8 

TSI-243-15  AL 

Tsinghai  Prov. 

0.  1 

-0.  5 

2.  9 

CRS-244-  14  AL 

Central  Russia 

2.  0 

1.  1 

2.  0 

PAK-  162-11  AL 

South  Pakistan 

1.  9 

4.  8 

4.  1 

PAK-179-06QC 

South  Pakistan 

3.  1 

4.  4 

4.  6 

PAK-179-10AL 

South  Pakistan 

4.  3 

5.  0 

3.  9 

*  =  Test  Area 
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•  TWN  -  SNNR  improvements  were  good  on  all  three  compon¬ 
ents.  The  Rayleigh  wave  improvements  were  better  than  the 
Love  wave  improvements. 

•  EKZ  -  SNNR  improvements  were  good  on  Love  waves  and 
excellent  on  Rayleigh  waves. 

In  general,  chirp  filter  Rayleigh  wave  SNNR  improvements 
were  better  than  Love  wave  improvements. 

The  standard  deviation  values  for  each  of  the  mean  SNNR  im¬ 
provements  indicate  there  is  considerable  variation  in  the  SNNR  improvements 
yielded  by  chirp  matched  filters.  For  example,  the  two  events  KAM-173-OOQC 
and  KAM- 180- 14 AL  have  SNNR  improvements  as  measured  on  LQ-T  differing 
by  more  than  4  dB.  A  possible  explanation  for  this  is  that  the  two  events  had 
different  source  mechanisms. 

The  chirp  length  data  is  summarized  in  Figures  V-3,  V-4  and 
V-5  for  the  transverse,  vertical,  and  radial  components  respectively.  The 
chirp  lengths  plotted  are  the  chirp  lengths  giving  the  best  improvements.  The 
corresponding  distances  are  great  circle  distances  in  kilometers  between 
event  epicenter  and  the  ALFA  array.  The  chirp  lengths  applied  to  presumed 
explosions  are  indicated  by  open  circles  while  those  applied  to  earthquakes 
are  indicated  by  dots.  A  least-mean- square -error  fit  was  made  for  the  data 
points  of  each  plot.  The  dotted  lines  100  seconds  above  and  100  seconds  be¬ 
low  this  line  represent  the  range  in  chirp  lengths  which  would  be  used  if  five 
chirps  were  applied  to  an  event,  the  chirp  length  increment  was  50  seconds, 
and  the  center  chirp  length  was  picked  from  the  least-mean- square-error  fit. 

For  each  component,  70  percent  or  more  of  the  data  points  fall 
within  the  100  second  bounds.  This  implies  that  the  least-mean-square-error 
fit  can  be  used  to  obtain  a  good  first  estimate  of  the  optimum  chirp  length  to 
apply  to  a  given  event.  The  range  in  epicentral  distance  for  each  region  is 
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CHIRP  FILTER  LENGTH  VS.  DISTANCE  (TRANSVERSE  COMPONENT) 


shown  on  Figures  V-3,  V-4,  and  V-5  to  indicate  the  dominant  data  source  for 
each  such  range. 

An  estimate  of  the  effectiveness  of  chirp  filters  was  made 
(using  LR-V  data)  by  computing  the  percentage  of  SNNR  improvements  less 
than  0  dB,  between  0  and  2  dB,  etc.  These  percentages  are:  15%  of  the  im¬ 
provements  were  less  than  0  dB,  21%  were  between  0  and  2dB,  24%  were  be¬ 
tween  2  and  4  dB,  30%  were  between  4  and  6  dB,  9%  were  between  6  and  8  dB, 
and  1%  were  greater  than  8  dB.  Thus  it  is  possible  to  say  that  the  probability 
that  the  output  trace  of  a  chirp  matched  filter  will  show  a  positive  improve¬ 
ment  is  0.  85  and  the  probability  that  the  improvement  will  be  between  2  and 
8  dB  is  approximately  0.  6  .  Therefore,  for  better  than  half  of  the  test  events 
to  which  chirp  filters  are  applied,  the  SNNR  improvements  will  be  large 
enough  to  aid  detection  studies. 

D.  COMPARISON  OF  REFERENCE  WAVEFORM  AND  CHIRP  MATCHED 
FILTER  RESULTS 

Those  regions  containing  four  or  more  test  events  (three  in  the 
case  of  EKZ)  were  selected  to  compare  reference  waveform  and  chirp  matched 
filter  results  on  a  regional  basis.  The  mean  SNNR  improvements  and  corres¬ 
ponding  standard  deviations  are  presented  in  Table  V-4  for  these  two  types  of 
filtering.  For  each  region,  the  computations  were  made  on  identical  sets  of 
test  events.  The  criteria  for  stating  that  one  type  of  filtering  performed 
better  than  the  other  was  that  the  better  filter  showed  a  higher  mean  improve¬ 
ment  and  a  lower  standard  deviation  of  the  improvement.  No  judgement  was 
made  in  cases  where  one  filtering  type  had  the  higher  mean  and  the  other  had 
the  lower  standard  deviation. 

The  results  of  this  are  summarized  as  follows: 

•  KAM  -  Chirp  filters  performed  better  than  reference  wave¬ 

form  filters  on  the  vertical  and  radial  Rayleigh  waves. 


•  CA  -  Chirp  filters  performed  better  than  reference  waveform 
filters  on  all  three  components. 

•  EKZ  -  Chirp  filters  performed  better  than  reference  waveforms 
on  the  Love  waves.  On  the  Rayleigh  waves,  chirps  were  better 

on  the  vertical  and  reference  waveforms  were  better  on  the  radial. 
Therefore,  their  performance  on  Rayleigh  waves  was  about  the 
same. 

•  TWN  -  The  results  in  this  region  were,  according  to  the  judge¬ 
ment  criteria,  indeterminant.  However,  on  the  Rayleigh  wave 
mean  improvements  for  reference  waveforms  are  almost  double 
those  of  the  chirps,  and  so  the  reference  waveforms  are  consid¬ 
ered  to  be  better. 

•  CSP  -  The  chirp  filters  performed  better  than  the  reference 
waveforms  on  the  vertical  Rayleigh  wave. 

•  SIR  -  The  chirp  filters  performed  better  than  the  reference 
waveforms  on  all  three  components. 

The  overall  relative  performance  of  chirp  and  reference  wave¬ 
form  matched  filters  is  illustrated  by  Figure  V-b.  This  figure  indicates  that 
57%  of  the  test  events  showed  higher  chirp  SNNR  improvements  than  reference 
waveform  improvements. 

One  other  comparison  between  chirp  and  reference  waveform 
matched  filter  performance  is  that  the  probability  that  a  matched  filter  will 
yield  a  SNNR  improvement  of  between  2  and  8  dB  is  0.45  for  reference  wave¬ 
forms  and  0.  60  for  chirps.  Thus,  it  is  concluded  that  chirp  matched  filters 
perform  slightly  better  than  reference  waveform  matched  filters. 

As  presently  conducted,  matched  filtering  is  fulfilling  the  goal 
of  increasing  the  ALPA  detection  capability  but  not  the  goal  of  providing  stable 
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estimates  of  SNNR  improvement  to  permit  the  computation  of  the  surface  - 
wave  magnitude  of  an  event  detected  only  by  a  matched  filter.  Inspection  of 
Appendix  A  reveals  that  of  128  events  from  1972  not  detected  on  the  bandpass- 
filtered  beam,  26  (20%)  were  detected  by  either  a  reference  waveform  or  a 
chirp  matched  filter.  Considering  the  1972  Eurasian  events  and  the  1972 
Kurile  Islands  -  Kamchatka  events  separately,  it  was  found  that  inclusion  of 
the  matched-filter  detections  did  not  change  the  90%  detection  thresholds  but 
did  change  the  50%  detection  thresholds  by  0.  3  magnitude  units  for  Eurasian 
data  and  0.4  magnitude  units  for  the  Kurile  Islands  -  Kamchatka  data. 

As  previously  described  in  this  section,  the  standard  deviations 
of  the  SNNR  improvements  are  so  large  that  the  use  of  an  "average"  SNNR 
improvement  to  compute  a  surface -wave  magnitude  is  not  very  meaningful. 
Therefore,  techniques  for  making  such  magnitude  estimates  from  matched 
filter  outputs  should  be  revised,  since  detections  made  on  matched  filter  out¬ 
puts  occur  only  when  the  match  is  good,  it  appears  that  using  a  high-side 
estimate  of  the  average  improvement  (e.  g.  ,  twice  the  mean  or  the  mean  plus 
one  standard  deviation)  would  give  a  better  surface-wave  magnitude  value. 

For  the  representative  case  of  3  dB  average  SNNR  improvement  and  3  dB 
variance,  this  would  mean  that  a  factor  of  2  should  be  used  to  get  the  surface- 
wave  magnitude  estimate  instead  of  a  factor  of  1.4,  which  would  lower  the 
surface -wave  magnitude  estimate  by  0.  I  units. 
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SECTION  VI 

S-WAVE  PROCESSING  RESULTS 


A.  INTRODUCTION 

S-wave  processing  was  performed  on  the  Kurile  Islands  - 
Kamchatka  earthquake  population  of  the  1972  data  base  to  resolve  an  apparent 
contradiction  appearing  in  the  report  of  the  preceeding  year  (Heiting,  et  al, 
1972).  In  that  report,  it  was  stated  that  the  90%  probability  of  S-wave  detec¬ 
tion  for  the  Kurile  Islands  -  Kamchatka  region  occurred  above  m  =5.0  How- 

b 

ever,  in  the  figure  accompanying  that  report,  this  90%  detection  level  appears 

to  be  at  about  m,  =  4.  4. 

b 

To  resolve  this  problem,  S-wave  processing  was  performed 
on  the  larger  population  of  events  available  from  the  1972  data  base.  Long- 
period  S-wave  beams  were  formed  for  103  events  having  epicenters  in  the 
Kurile  Islands  and  Kamchatka  regions,  using  the  apparent  horizontal  S-wave 
velocity  appropriate  to  each  epicentral  distance.  The  values  of  S-wave  velo¬ 
city  were  taken  from  a  plot  of  S-wave  apparent  horizontal  velocity  as  a  func¬ 
tion  of  epicentral  distance  (Texas  Instruments  Incorporated,  1964).  Band¬ 
pass  filtered  (0.  025  to  0.  055  Hz)  beams  were  formed  for  the  rotated  trans¬ 
verse,  vertical,  and  radial  components.  The  detection  criteria  were  that 
the  signal  be  at  least  3  dB  above  the  noise  and  that  the  signal  arrived  within 
+.10  seconds  of  the  predicted  arrival  time.  The  amplitude /period  data  were 
measured  on  the  bandpass-filtered  component  which  showed  the  largest  S- 
wave  amplitude. 


B.  LONG-PERIOD  S-WAVE  DETECTION  THRESHOLD  ESTIMATE  FOR 
ALPA 

The  histogram  in  the  upper  portion  of  Figure  VI- 1  shows  the 
total  number  of  events  processed  at  each  body-wave  magnitude  for  the  Kurile 
Islands  -  Kamchatka  earthquake  populations.  These  are  subdivided  (shaded 
and  unshaded  portions  of  the  histogram)  into  summer  and  winter  populations. 
The  lower  graph  shows: 

•  detection  percentages  as  a  function  of  body-wave  magnitude 
for  the  total  population, 

•  detection  percentages  for  the  summer  population, 

•  detection  percentages  for  the  winter  population. 

The  data  clearly  indicates  that  the  90%  probability  of  S-wave 

detection  occurs  just  below  m  =  4.  5.  This  also  holds  t.-ue  for  the  summer 

b 

and  winter  populations.  However,  it  is  worthy  of  note  that  a  much  higher 

percentage  of  summer  events  (30%)  with  m,  values  below  4.  5  were  detected 

b 

than  of  winter  events  (7%).  This  may  be  attributed  to  the  lower  RMS  noise 
level  of  the  summer  months. 

C.  DISCRIMINATION  BY  MEANS  OF  LONG-PERIOD  S- WAVES 

The  S-wave  amplitude /period  data  for  the  36  detected  events 
are  displayed  in  Figure  VI-2.  For  comparison,  the  corresponding  data  for 
the  few  presumed  explosions  which  had  detected  S-waves  are  indicated  by 
solid  triangles.  The  values  of  A/ T  were  normalized  to  a  body-wave  magnitude 
of  5.  0  and  an  epicentral  distance  of  20°,  following  the  general  method  used  by 
Evernden  (1969).  The  normalization  computation  may  be  summarized  as: 

A/T 

Norm. 

(Heiting,  et  al,  1972) 


'io0- 
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m  +  4.  74 
b 


A/T 
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The  normalized  data  show  complete  separation  by  a  factor  of 
two  between  the  earthquake  and  presumed  explosion  populations.  Also,  there 
is  a  factor  of  ten  separation  between  the  average  A/T  values  of  the  two  pop¬ 
ulations.  Therefore,  S-wave  A/T  normalized  values  are  a  good  earthquake- 
presumed  explosion  discriminant  for  m  =  4.  5  or  greater. 
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SECTION  VII 

ALPA  SURFACE  WAVE  DETECTION  CAPABILITY 


Estimates  of  ALPA  surface  wave  detection  capability  were  de¬ 
termined  for  each  region  described  in  Section  IV  with  the  exception  of  the 
Taiwan  region,  which  did  not  have  a  large  enough  event  population.  In  each 
case,  this  estimate  was  obtained  by  plotting  the  percentage  of  earthquakes 
for  which  surface  waves  were  detected  as  a  function  of  body-wave  magnitude. 
The  histogram  in  the  upper  portion  of  each  of  the  Figures  V1I-1  through  V1I-7 
describes  the  portion  of  the  1970-1972  data  base  belonging  to  the  given  region. 
These  histograms  were  used  to  compute  the  incremental  detection  probabili¬ 
ties  which  appear  in  the  lower  portion  of  each  of  the  figures.  The  events  used 
are  designated  D  for  detected  and  ND  for  not  detected  in  Appendix  A. 

The  criteria  for  determining  whether  detection  was  achieved 
for  an  event  are: 


•  A  peak  in  any  output  trace  3  dB  above  any  other  peak  in  a  20 
minute  time  gate  centered  at  the  expected  peak  occurrence 
time. 

•  A  peak  which  occurs  within  ±  180  seconds  of  the  expected  peak 
occurrence  time. 

It  should  be  noted  that  these  criteria  are  not  absolute.  Occasionally,  an  event 
was  listed  as  detected  which  did  not  meet  both  of  these  criteria.  Peaks  could 
be  less  than  3  dB  above  other  peaks  in  the  20  minute  time  gate  and  still  be 
identified  as  signals  from  their  dispersion  characteristics.  Also,  some  peaks, 
from  Central  Asian  events  in  particular,  appeared  later  than  the  second  cri¬ 
terion  allows,  bat  were  still  recognized  as  signals.  Although  we  do  not  have 
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FIGURE  VII- 1 

SURFACE  WAVE  DETECTION  DATA  FOR  KAMCHATKA  REGION 
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SURFACE  WAVE  DETECTION  DATA  FOR  CASPIAN  SEA  REGION 
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a  quantitative  figure,  we  believe  that  the  false  alarm  rate  associated  with  this 
detection  scheme  is  very  low  (<  1  %). 


By  region,  the  results  of  this  analysis  are: 

Kamchatka  Region:  The  90%  detection  level  is  at  m  =  4.  1  . 

b 

The  notch  in  the  detection  curve  at  m  =4.4  (indicated 

b 

by  an  arrow)  was  not  considered  significant,  since  there 

are  only  two  events  at  this  m  increment,  while  there 

b 

are  seven  at  m  =  4.  3  and  5  at  m  =  4.  2,  all  of  which 
D  b 

were  detected. 


Kurile  Islands  Region:  The  90%  detection  level  is  between 

m  =  4.  1  and  m,  =  4.  2. 
b  b 

Central  Asia  Region:  The  90%  detection  level  is  just  below 

m  =  4.  4. 
b 


Caspian  Sea  Region:  The  90%  detection  level  is  at  m  =  4.  5 

b 

Southern  Iran  Region:  The  90%  detection  level  is  at  m  =  4.  3 

b 

Greece-Turkey  Region:  The  90%  detection  level  is  at  m  =  4.  5. 

b 

The  notch  in  the  detection  curve  at  m^  =4.8  (indicated  by 
an  arrow)  was  not  considered  significant,  since  it  was 
caused  by  an  event  which  was  reported  by  only  three  sta¬ 


tions. 

fore, 


U1  tnese,  only  one  reported  a  value  of  m  .  There - 

b 

the  mb  of  this  event  was  not  considered  to  be  reliabl 


Eastern  Kazakh  Test  Area:  The  90%  detection  level  for  pre¬ 
sumed  explosions  from  this  area  is  just  below  m  =  5.  5. 

b 

Even  though  the  event  population  for  this  area  is  small, 
it  was  felt  that  this  detection  level  estimate  could  be  made, 
since  all  members  of  the  population  were  presumed  ex¬ 
plosions  and  all  were  from  one  small  area. 


These  results  show  that,  insofar  as  the  scope  of  the  data  base 
and  the  90%  detection  level  are  concerned,  only  four  general  areas  need  be 
considered.  These  are  the  Kamchatka  -  Kurile  Islands  region,  the  Central 
Asia  region  including  Southern  Iran,  the  region  encompassing  Southeastern 
Europe  and  the  Caspian  Sea  area,  and  the  Eastern  Kazakh  test  area. 

In  addition  to  the  regional  detection  levels  described  above, 
histograms  and  their  corresponding  detection  curves  were  constructed  for 
winter  (January  through  March)  and  summer  (June  through  August)  suites  of 
events.  These  are  presented  in  Figures  VII-8  and  VII-9.  The  90%  detection 
level  for  the  winter  suite  lies  at  mb=  4.  5,  while  the  90%  detection  level  for  the 
summer  suite  lies  at  =  4.4.  This  difference  of  0.  1  magnitude  unit  in  the 
detection  levels  can  be  ascribed  to  the  bursts  of  long-period  non-propagating 
noise  occuring  during  the  winter  months  (Section  III),  which  would  tend  to 
obscure  events  which  would  otherwise  have  been  detected. 

A  comparison  of  the  results  yielded  by  this  analysis  with  those 
reported  for  the  preceeding  year  (Heiting,  et  al,  1972)  shows  that  the  90%  de¬ 
tection  level  is  lower  this  year  by  0.  1  m  units  for  both  the  Central  Asia  and 

b 

Kurile  Islands  -  Kamchatka  regions.  This  is  probably  due  to  the  greater 
precision  possible  this  year  because  of  a  larger  data  base,  and  to  the  exclu¬ 
sion  this  year  of  events  from  the  Central  Asia  region  which  more  properly 
belong  in  the  Southeastern  Europe  or  Caspian  Sea  regions. 
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SURFACE  WAVE  DETECTION  DATA  FOR  SUMMER  EVENT  SUITE 
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SECTION  VIII 


BEHAVIOR  OF  STANDARD  DISCRIMINANTS 

A.  MEASUREMENT  OF  M 

s 

Surface  wave  magnitudes  (MJ  were  computed  for  all  of  the 

events  which  were  detected  and  listed  in  Appendix  A.  Whenever  possible, 

values  of  M  were  computed  for  all  three  components.  For  events  which 
s 

were  not  detected,  an  upper  bound  for  the  corresponding  surface  wave  mag¬ 
nitude  was  computed  from  the  largest  peak-to-peak  noise  amplitude  occur¬ 
ring  within  the  signal  gate.  These  upper  bounds  are  designated  by  the 
symbol  )  preceeding  the  M  value  in  Appendix  A.  The  surface  wave  mag¬ 
nitude  was  determined  for  each  event  by  the  formula: 

M^  =  log  A/T  +  1.  66  log  A 

where: 

A  is  the  largest  peak-to-peak  amplitude  in  millimicrons 
T  is  the  period  in  seconds  in  the  neighborhood  of  the  peak 
A  is  the  epicentral  distance  in  degrees. 

All  measurements  were  made  on  beamsteered  traces  which 
had  been  bandpass  filtered  (0.  025  -  0.  055  Hz  passband).  The  maximum  peak 
amplitudes  usually  occurred  near  25  second  periods.  However,  for  some  of 
the  more  distant  events,  the  maximum  peak  amplitudes  were  measured  at 
periods  between  28  and  32  seconds.  To  estimate  the  effect  of  the  period  at 
which  the  amplitude  was  measured  on  the  resulting  value  of  M  ,  A/T  values 
for  58  events  were  measured  twice;  once  in  the  period  range  20-25  seconds 
and  once  in  the  period  range  26-32  seconds.  The  results  indicated  that  the 
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values  of  for  Love  waves  were  essentially  the  same  for  the  two  period 

ranges,  while  the  values  of  for  Rayleigh  waves  as  measured  in  the  period 

range  26-32  seconds  may  be  as  much  as  0.  1  M  units  lower  than  those  mea- 

s 

sured  in  the  first  period  range.  Since  the  accuracy  of  measurement  only 

approaches  0.  1  M  units,  all  M  values  reported  were  measured  on  the  lar- 
s  s 

gest  amplitude  present  on  the  0.  025  to  0.  055  Hz  bandpass  filtered  trace  with¬ 
out  further  concern  for  the  period  at  which  the  amplitude  was  measured. 


B.  M  -  m,  AS  A  DISCRIMINANT 
s  b 

The  M^-  m^  data  for  each  region  are  presented  in  Figures 
VIII- 1  through  VIII-4.  Only  three  regions  had  large  enough  event  populations 
to  allow  comparison  with  the  presumed  explosion  population.  One  other  re¬ 
gion,  Taiwan,  is  presented;  even  though  the  event  population  is  small,  all  of  it 

events  were  detected.  Each  figure  consists  of  two  plots  -  one  for  M  mea- 

s 

sured  on  Love  waves  and  one  for  measured  on  Rayleigh  waves.  In  each 
plot,  the  corresponding  data  for  presumed  explosions  are  presented  for  pur¬ 
poses  of  comparison.  On  these  plots,  a  circle  represents  an  earthquake  M  - 

s 

m^  data  point,  and  a  triangle  represents  a  presumed  explosion  M  -  m  data 

s  b 

point.  A  vertical  line  below  one  of  these  symbols  indicates  an  upper  bound 

on  M  for  a  non-detected  event, 
s 


A  least-mean-square-error  fit  was  made  for  each  data  set.  To 

avoid  bias  at  lower  values  of  m^  where  only  those  events  having  a  relatively 

high  M  were  detected,  the  fit  was  made  only  on  those  values  of  m  which 
s  b 

were  at  or  above  the  100  per  cent  detection  level.  Data  points  at  values  of  m 

b 

greater  than  5.  5  were  also  excluded  from  the  fit,  since  above  this  m  value 

b 

the  fit  assumes  a  steeper  slope  (Tsai,  1972).  These  least-square  fits  are 
shown  as  solid  lines  over  their  intervals  of  definition.  The  dashed-line  ex¬ 
tension  of  these  fits  are  presented  as  an  aid  in  comparing  the  earthquake  pop¬ 
ulation  fit  with  the  presumed  explosion  fit.  With  the  exception  of  the  TWN 
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region  data  and  the  KUR  region  LR-V  data  (which  have  slopes  less  than  1.  0), 
the  slopes  of  the  earthquake  populations  are  between  1.  1  and  1.4  and  are 
therefore  parallel  or  near-parallel  to  the  LQ-T  and  LR-V  fits  of  the  pre¬ 
sumed  explosion  data. 

In  general,  the  plots  show  good  separation  between  earthquake 
and  presumed  explosion  populations.  One  presumed  explosion,  EKZ-345-04AL, 
has  an  abnormally  high  value  of  for  both  Love  and  Rayleigh  waves.  (This 

event  is  indicated  by  #  on  the  plots.)  However,  the  SD  AC/ LAS  A  weekly 
summary  reports  this  presumed  explosion  as  a  presumed  double  explosion 
with  8  second  separation.  If  this  is  the  case,  the  high  M  values  shown  by 
this  event  can  be  explained  by  constructive  interference  of  the  two  events. 

The  poor  separation  shown  by  the  Kurile  Islands/  Kamchatka  events  as  report¬ 
ed  last  year  (Heiting,  el  al,  1972)  now  appears  to  be  due  to  the  events  from 
the  southern  Kurile  Islands.  The  Kamchatka  region  events  show  complete 
separation  from  the  presumed  explosion  population  (with  the  exception  of  the 
presumed  explosion  EKZ-345-04AL,  as  described  above).  The  Central  Asia 
and  Taiwan  event  populations  also  show  complete  separation  from  the  pre¬ 
sumed  explosion  population.  Comparison  of  the  least-mean-square-error 

fits  indicates  that,  in  every  region,  Love  wave  M  -  m,  is  a  better  discrimi- 

s  b 

nant  than  Rayleigh  wave  M  -  m  . 

s  b 

Using  only  those  earthquakes  having  values  at  or  above  the 
100  per  cent  detection  level  and  below  5.  6,  M  -  m  relationships  were  com- 
puted  for  the  entire  body  of  data.  By  also  computing  the  95  per  cent  confidence 
limits  for  these  relationships,  it  was  possible  to  estimate  the  error  involved 
in  the  computation.  The  resulting  M  -  m  relationships  are: 


for  Love  wave  M  : 

s 

for  Rayleigh  wave  M  : 


=  (l.Z  +  0.  2)m^  -  (2,  1  i  1,  1)  (163  events) 
-  (1*  2+0.  2)m^  -  (2.  2  +  0,  9)  (205  events) 
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These  relationships  are  clearly  different  than  the  Gutenber g- Ri chte r  relation¬ 
ship: 

M  =  1.  59  m  -  3.  97 
s  b 

and  strongly  suggest  that  the  M  -  m  slope  in  the  region  5.  6  4  m  4  4.2  is 

s  b  b 

significantly  lower  than  that  observed  at  higher  magnitudes. 

C.  AD  AND  AR  AS  DISCRIMINANTS 

The  parameter  AR  is  related  to  the  total  Rayleigh  wave  energy 
in  a  seismic  event.  It  was  introduced  by  Brune,  Espinosa,  and  Oliver  (1963). 
It  has  been  used  by  Evernden  (1969),  who  also  studied  AE,  the  corresponding 
parameter  for  Dove  wave  energy.  The  AR  and  AD  parameters  as  used  here 
were  computed  by  summing  the  absolute  values  (in  millimicrons)  of  the  data 
points  within  the  signal  gate  beginning  at  the  expected  arrival  time  of  the  sig¬ 
nal  and  extending  throughout  the  expected  time  of  duration  of  the  signal.  The 
results  were  scaled  as  described  earlier  (Heiting,  et  al,  197Z). 

Only  two  regions,  Kamchatka  and  Central  Asia,  contained 
enough  data  to  make  analysis  of  AD  and  AR  meaningful.  Figures  Vlll-5  and 
VIII-6  present  these  data.  Open  circles  represent  values  of  AD;  solid  circles 
represent  values  of  AR.  For  comparison,  presumed  explosion  values  of  AD 
and  AR  are  plotted,  using  open  and  solid  triangles  respectively. 

Examination  of  the  data  reveals  that  the  parameters  AD  and 
AR  have  a  fair  discrimination  capability.  While  there  is  some  overlap  of  the 
earthquake  and  presumed  explosion  populations,  the  average  separation  is  on 
the  order  of  a  factor  of  ten.  The  highest  AD  and  AR  values  for  a  presumed 
explosion  were  measured  on  EKZ  -  345 -04  AD.  If  this  event  is  removed  from 
the  presumed  explosion  data,  for  reasons  stated  earlier,  the  Central  Asia 
AD  and  AR  values  will  show  no  overlap  with  the  presumed  explosion  values. 
The  overlap  of  the  Kamchatka-presumed  explosion  AR  data  will  be  greatly 
reduced  and  there  will  be  no  overlap  for  the  corresponding  AD  data.  In  detail, 
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the  separation  between  the  lowest  earthquake  value  and  the  highest  presumed 
explosion  value  will  be  a  factor  of  1.8  for  Kamchatka  AL,  3.  5  for  Central 
Asia  AL.,  and  1.5  for  Central  Asia  AR.  Therefore,  AL  appears  to  be  a 
better  discriminant  than  AR. 


SECTION  IX 
CONCLUSIONS 

A.  MAJOR  RESULTS 

Summarized  below  are  trie  major  results  of  each  of  the  areas 

of  evaluation: 

1.  Noise  Analysis 

•  There  appears  to  be  a  fairly  constant  background  RMS 
noise  level  throughout  the  year  at  ALPA  in  the  range  of 
7-10  m/i  . 

•  For  the  most  part,  the  higher  RMS  noise  levels  (>14  my) 
appear  to  be  due  to  sudden  increases  in  long-period  non¬ 
propagating  noise  superimposed  on  the  background  noise 
level.  Only  a  few  of  the  high  RMS  noise  levels  may  be 
due  to  storm-generated  noise. 

•  The  source  azimuths  of  microseismic  noise  as  recorded 
at  ALPA  rarely  coincide  with  azimuths  to  the  area  of  in¬ 
terest. 

2.  Matched  Filter  Studies 

•  There  is  considerable  variation  in  SNNR  improvement 
within  a  given  region  for  both  reference  wavefo  m  and 
chirp  matched  filter 

•  The  standard  deviations  are  so  large  that  mean  SNNR 
improvement  values  are  almost  meaningless. 
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•  Chirp  matched  filters  appear  to  be  slightly  more  effec¬ 
tive  than  reference  waveform  matched  filters. 

•  Although  the  use  of  matched  filters  does  not  appreciably 
affect  the  90%  detection  thresholds,  it  increases  the  num 
ber  of  events  detected  and  lowers  the  50%  detection 
thresholds. 


S-Wave  Processing  Results 


4. 


•  The  90%  S-wave  detection  threshold  is  just  below  m  = 

4.  5  for  the  Kamchatka  -  Kurile  Island  events. 

•  A  higher  percentage  of  summer  events  (30%)  with  m 
below  this  90%  detection  threshold  was  detected  than  of 
winter  events  (7%). 

•  F°r  events  Wlth  mb  =  4.  5  or  greater,  S-wave  A/ T  values 
are  a  good  discriminant. 


■^^PA  Surface  Wave  .Detection  Capability 


By  region,  the  90%  detection  level  (with  a  corresponding 
low  «1%)  false  alarm  level)  for  surface  waves  occurs  at 


Kamchatka 
Kurile  Islands 
Central  Asia 
Caspian  Sea 
Southern  Iran 
Greece -Turkey 
Eastern  Kazakh 


between  m  =  4.  1  and  m  =4.2 
D  b 

between  m  =  4.  3  and  m,  =  4.  4 
D  b 

mb  =  4.5 

m  =4.3 
b 

m  =4.5 
b 

mb  =  5,5  (Presumed  Explosions  Only) 


The  90%  detection  levels  for  the  winter  and  summer  event 
suites  are  at  mb  =  4.5  for  the  winter  and  at  m  =  4.  4  for 
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the  summer.  This  difference  in  detection  levels  is  be¬ 
lieved  to  be  due  to  the  slightly  higher  noise  levels  of  the 
winter  months. 

5.  Behavior  of  the  Standard  Discriminants 

•  As  long  as  the  peak-to-peak  amplitude  of  a  surface  wave 
was  measured  at  the  highest  amplitude,  the  resulting 
value  of  did  not  appear  to  depend  upon  the  period  of 
the  waveform  at  which  this  amplitude  was  measured. 

•  With  the  exception  of  one  presumed  explosion,  which  may 

actually  have  been  two  explosions,  the  M  -  m.  discrimin- 

s  b 

ation  method  achieved  complete  separation  of  earthquakes 

and  presumed  explosions.  The  Love  wave  M  -  m  appears 

s  b 

to  be  a  better  discriminant  than  the  Rayleigh  wave  M  -  m  . 

s  b 

•  Using  all  available  earthquake  data  with  mb  values  equal 
to  or  greater  than  the  100%  detection  level  and  less  than 
5.  6,  the  following  M  -  m^  relationships  were  derived: 

for  Love  wave  M  :  M  =  (1.2+0.  2  )m  -(2.  1  +  1.1) 

s  s  —  b 

for  Rayleigh  wave  M  :  M  =  (1. 2  ±  0.  2)m  -  (2.  2  +  0.  9) 

S  S  D 

The  error  estimates  for  the  slope  and  intercept  were 
measured  from  the  95%  confidence  limits  of  the  computed 
least-mean- square -error  fits  to  the  data. 

•  AL  and  \R  were  not  as  successful  in  discrimination  as 

M  -  m.  .  With  the  removal  from  the  data  set  of  the  pre- 
s  b  K 

sumed  double  explosion,  separation  between  the  earthquake 
and  presumed  explosion  populations  was  complete  except 
for  the  Kamchatka  area  AR  data.  AL  data  showed  better 
separation  than  AR  data. 
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B. 


SUGGESTIONS  FOR  FUTURE  ANALYSIS 


The  following  areas  should  be  investigated  in  any  future  ana¬ 
lysis  of  ALP  A. 

•  The  indicated  upward  trend  of  the  RMS  noise  level  from  day 
241  to  day  361  of  1972  should  be  investigated  to  more  fully 
determine  if  it  is  different  from  the  corresponding  time  period 
of  1971.  This  would  require  more  noise  samples  taken  in  the 
above  mentioned  period  and  an  extension  of  this  period  into 
1973. 

•  More  data  for  reference  waveform  matched  filters  and  chirp 
matched  filters  should  be  compiled  to  make  possible  a  more 
thorough  investigation  of  their  regional  characteristics. 

•  More  events  from  Central  Asia  should  be  processed  to  make 
possible  the  subdivision  of  this  region.  In  particular,  the 
differences  in  detection  of  Western  Sinkiang  events  and  Hindu 
Kush  events  should  be  investigated.  At  present,  the  small 
number  of  events  in  these  areas  makes  it  impossible  to  do 
more  than  note  that  most  events  from  Western  Sinkiang  were 
detected,  but  few  from  the  Hindu  Kush  area  were  detected. 


C.  GENERAL  CONCLUSIONS  OF  THE  THREE -YEAR  ALPA  ARRAY 
EVALUATION 

With  the  termination  of  the  three -year  evaluation  of  the  detec¬ 
tion  and  discrimination  capabilities  of  ALPA,  the  following  conclusions  can 
be  made  about  its  characteristics  and  performance: 

•  Signal  similarity  across  the  array  generally  is  good.  As  ex¬ 

pected,  similarity  across  the  full  19-eiement  array  is  less 
than  that  .'.cross  the  limited  9-element  array.  The  a/erage 
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signal  correlation  coefficient  for  the  vertical  component  is 
0.  84  for  the  full  array  and  0.  93  for  the  limited  array. 

The  beamsteer  signal  attenuation  averages  about  2  dB  across 
the  full  array  and  about  1  dB  for  the  seven-element  hexagonal 
subarray  on  all  three  components. 

The  noise  field  at  ALPA  is  characterized  by  a  fairly  constant 
background  RMS  noise  level  of  7  to  1 0  m/i  (on  a  single  channel) 
which  is  punctuated  during  the  winter  months  by  bursts  of  long- 
period  non-propagating  noise.  These  bursts  can  temporarily 
double  or  triple  the  RMS  noise  level  and  hence  decrease  detec¬ 
tion  capability.  Propagating  storm-generated  noise  is  confined 
to  a  narrow  band  around  18  seconds  and  variations  in  this  peak 
do  not  significantly  affect  detection  capability.  The  source  azi¬ 
muths  of  propagating  microseismic  noise  rarely  coincide  with 
azimuths  to  the  area  of  interest.  Noise  levels  are  essentially 
the  same  on  all  three  components. 

Noise  reduction  achieved  by  beamsteering  is  very  close  to  Mn; 

hence  output  beam  RMS  noise  levels  usually  are  between  1.  5 
and  2.  5  m^i. 

The  ALPA  noise  is  not  time  stationary;  substantial  variations 
in  wave  number  structure  have  been  observed  at  the  micro- 
seismic  peak  in  a  two-hour  period.  Unless  the  design  noise  is 
within  a  very  few  hours  of  the  data  to  which  a  multichannel 
filter  is  to  be  applied,  there  is  no  advantage  of  a  multichannel 
filter  over  beamsteering. 

The  matched  filter  studies  indicate  that,  in  general,  chirp 
matched  filters  perform  slightly  better  than  reference  waveform 
matched  filters  in  that  they  yield  essentially  the  same  mean 
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SNNR  improvements  as  the  reference  waveforms  but  have  less 
variation  in  improvement  among  test  events.  Since  matched 
filters  decreased  the  number  of  otherwise  undetected  events 
by  about  20%,  thereby  lowering  the  50%  detection  levels,  they 
are  of  value  in  event  detection  studies. 

•  Two -component  beamforming  as  performed  in  1971  yielded 
SNNR  gains  over  the  one  component  beam  of  only  one  to  two 
dB  in  the  bandpassed  output  beam.  Therefore,  two -component 
beamforming  was  considered  unsatisfactory  as  a  signal  en¬ 
hancement  technique  and  was  not  used  in  the  1972  evaluation. 

•  The  S-wave  90%  detection  threshold  was  determined  to  be  at 

m.  =4.5  for  Kurile  Islands  -  Kamchatka  events  and  at  m,  = 
b  b 

5.  5  for  Central  Asian  events.  The  S-wave  A/ T  discriminant 
is  good  for  events  having  m^  values  at  or  above  these  detection 
thre  sholds. 


The  surface  wave  90%  detection  levels  were  determined  to  be 

at  m,  =  4.  1  for  Kurile  Islands  -  Kamchatka  events,  at  m.  = 
b  b 

4.  4  +_  0.  1  for  Eurasian  events,  and  at  m^  =  5,5  for  Eastern 

Kazakh  presumed  explosions.  The  winter  suite  90%  detection 

level  was  found  to  be  at  rn  =4.5  and  the  summer  90%  detec- 

b 

tion  level  at  m,  =  4.4.  This  difference  is  believed  to  be  due 
b 

to  the  higher  RMS  noise  levels  occurring  during  the  winter 
months.  Detection  levels  for  a  nine -element  subarray  are 
only  0.  1  to  0.  2  magnitude  units  above  those  for  the  full  array. 


The  best  earthquake -pre sumed  explosion  discriminant  appears 

to  be  the  M  -  m,  relationship.  The  M  -  m,  relationship  de- 
s  b  s  b 

termined  from  Love  wave  energy  is  a  slightly  better  discrim¬ 
inant  than  the  corresponding  relationship  determined  from  Ray¬ 
leigh  wave  energy.  The  AL-m^  and  AR-m^  discriminants, 
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while  inferior  to  the  M  -m,  discriminants,  are  useful  in  eartn- 

s  b 

quake  presumed  explosion  discrimination  studies.  AL-m,  is 

b 

a  better  discriminant  than  AR-m  . 

b 

Using  only  those  earthquakes  having  m,  values  at  or  above  the 

b 

100%  detection  level  and  below  m,  =  5.  6,  the  following  M  -m 

b  s  b 

relationships  were  determined: 


for  Love  wave  M  :  M  =  (1.2  +  0.  2) m.  -  (2.1  +.1.1) 

s  s  b 

for  Rayleigh  wave  M  :  M  =  (1.2  +  0.  2)m.  -  (2.2  +0.9) 

The  error  estimates  for  the  slope  and  intercept  were  deter¬ 
mined  from  the  95%  confidence  limits.  These  relationships 
are  clearly  different  than  the  Gutenberg -Richter  relationship 

M  =  1.  59  m.  -  3.  97  and  indicate  that  the  M  -  m,  slope  in  the 
s  b  s  b 

4.  2  <  <5.6  range  is  lower  than  that  at  higher  magnitudes. 
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APPENDIX  A 

THE  COMBINED  197  0  -  1972  DATA  BASE 


All  events  used  in  this  evaluation  of  ALP  A  are  listed  on  the 
following  pages.  The  events  are  listed  in  chronological  order  starting  on 
May  17,  1970  and  ending  on  December  10,  1972.  The  parameters  listed 
are:  event  name,  month  of  occurrence  (Mo.  ),  day  of  occurrence  (Day), 
origin  time,  latitude  (Lat/°N),  longitude  (Long/°E),  body  wave  magnitude 
(m  )3  Rayleigh  surface  wave  magnitude  ((MS)R)»  Eove  surface  wave  mag¬ 
nitude  ((M  )  ),  region  as  defined  in  Section  IV  (Reg),  depth  (D),  detection 

(Det.  ),  number  of  sites  used  in  processing  (NS),  and  information  source  (IS). 


follows : 


The  symbols  accompanying  some  of  these  parameters  are  as 

preceeding  m^  value  —  m^  has  been  recalculated  using  tele- 

seismic  data  only. 


preceeding  M  value  -  upper  bound  of  Mg  for  a  non-detected 

event. 


M 


preceeding  M 


s 


value 


detection  by  reference  waveform  match 
ed  filter. 


C  preceeding  K4  value  -  detection  by  chirp  matched  filter. 

X  preceeding  information  source  -  presumed  explosion 

in  any  column  -  no  information  available. 


A-  1 


u  -i  iv  rmiV'iif-m-rr 


The  information  source  code  is: 

P  -  parameters  taken  from  Preliminary  Determination  of  Epi¬ 
centers  Monthly  Summary 

I  -  parameters  taken  from  International  Seismological  Month 
verified  event  list 

J  -  parameters  taken  from  International  Seismological  Month 
unverified  event  list 

L  -  parameters  taken  from  SDAC/LASA  Weekly  Summary 

N  -  parameters  taken  from  NORSAR  Seismic  Event  Summary. 
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COMPLETE  LIST  OF  EVENTS  USED  IN  ALPA  EVALUATION 

(PAGE  1  OF  20) 
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APPENDIX  B 

DISPOSITION  OF  EVENTS  PROPOSED  FOR  ROUTINE  PROCESSING 


A  breakdown  of  the  disposition  of  events  is  presented  in  the 

table  below. 


Month 

Total  Number  of 
Events  Proposed 
For  Processing 

Processing 

Completed 

Not  Completed  Due 
To  Interference 

Not  Completed  Due 
To  Tape  and  Data 
Problems 

January 

81 

67% 

17% 

16% 

February 

131 

57% 

34% 

9% 

March 

84 

48% 

18% 

34% 

June 

122 

65% 

19% 

16% 

July 

112 

62% 

12% 

26% 

August* 

146 

44% 

27% 

29% 

*  The  processing  of  August  events  was  not  completed  due  to  time  re¬ 
straints.  The  percentages  are  based  on  the  number  of  events  for 
which  processing  was  attempted  (132). 

The  term  "interference"  as  used  here  is  defined  in  Quarterly  Report 
No.  3. 

Tape  and  data  problems  include  parity  errors,  uncorrectable  spikes, 
and  data  not  on  tape. 
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APPENDIX  C 


NOISE  ANALYSIS  PARAMETERS 
1972  NOISE  SAMPLES 


Start  Time 

No.  of  Sites 
Available 

Frequency 

Azimuth 

Velocity 
km/  sec 

1 

08.  20.  00 

11 

.  055 

135° 

3.  9 

12 

22.  20.  00 

8 

.  063 

225° 

3.  6 

21 

04.  50.  00 

13 

.  059 

148° 

3.  5 

31 

14.  00.  00 

11 

.  063 

237° 

3.  8 

41 

22  00.  00 

15 

.  063 

127° 

3.6 

51 

23.  15.  00 

15 

.  059 

178° 

3.  5 

61 

08.  00.  00 

10 

.  051 

19° 

4.  0 

71 

14.  40.  00 

13 

.  051 

145° 

3.  8 

81 

21.  20.  00 

11 

.  059 

172° 

3.  6 

90 

01.  30.  00 

9 

.  063 

189° 

3.7 

101 

19.20.  00 

13 

.  055 

27° 

3.  9 

111 

23.  00.  00 

10 

.  055 

211° 

3.  8 

121 

23.  15.  00 

10 

.  059 

140° 

3.  6 

131 

00.  55.  00 

9 

.  055 

141° 

3.  9 

142 

09.45.  00 

7 

.  063 

222° 

3.8 

152 

02.  20.  00 

10 

.  051 

207° 

3.  9 

1 6l 

22.  30.  00 

12 

.  059 

143° 

3.  8 

172 

20.  55.  00 

13 

.  051 

232° 

3.  6 

181 

02.  25.  00 

11 

.  055 

163° 

3.  9 

191 

00.  45.  00 

7 

.  055 

288° 

3.8 

201 

00.  15.  00 

15 

.  051 

141° 

3.  6 

211 

15.  00.  00 

16 

.  055 

140° 

3.  6 

220 

18.  38.  00 

9 

.  059 

209° 

3.  6 

232 

20.  40.  00 

11 

.  059 

288° 

3.  7 

241 

01.  10.  00 

16 

.  059 

215° 

3.6 

250 

10.  30.  00 

12 

.  055 

188° 

3.  5 

261 

08.  20.  00 

11 

.  059 

288° 

3.  8 

272 

14.  10.  00 

11 

.  047 

89° 

3.  7 

281 

10.  50.  00 

12 

.  059 

180° 

3.  5 

291 

20.  30.  00 

14 

.  063 

139° 

3.8 

C-l 


Day 


Start  Time 


No.  of  Sites 
Available 


Azimuth 


301 

22.40.  00 

11 

.  059 

192° 

311 

08.40.  00 

15 

.  055 

48° 

321 

10.  20.  00 

10 

.  059 

143° 

331 

20.  20.  00 

15 

.  063 

141° 

341 

17.  50.  00 

15 

.  055 

0° 

352 

19.45.  00 

9 

.  055 

68° 

361 

15.  30.  00 

16 

.  063 

135° 

